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We analyze the adiabatic magnetization of ferromagnetic clusters in an intermediate coupling 
regime, where the anisotropic potential is comparable to other energy scales. We find a non- 
monotonic behavior of the magnetic susceptibility as a function of coupling with a peak. Coriolis 
coupling effects are calculated for the first time; they reduce the susceptibility somewhat. 



I. INTRODUCTION 

Magnetic-properties of a wide variety of ferromagnetic clusters, Fe, Ni and Co in transition metals involving 3d 
electronsalffitl, and Ru, Rh and Pd associated with 4d electronstna, are studied using the Stern-Gerlach technique. 
The observed deflection profile, caused by the interaction of magnetization and the gradient of field strength, provides 
information about these intrinsic magnetic moment of the cluster as well as its coupling to the other degrees of freedom 
in clusters. The size of cluster is small enough to be regarded as a single domain system, and the electrons involved 
form a single giant magnetic moment; we will call it the super-electron spin in this paper. It is quite interesting to 
see that, in such a low dimensional system, some elements form ferromagnetic clusters, even though the bulk material 
is nonmagnetic. Also, the magnetization is strongly dependent on the number of atoms in the clusters. Therefore, 
it is important in the present stage of study to establish a method of analysis for extracting the intrinsic magnetic 
moment from the observed deflection profile, and this is the motivation for the present work. 

Let us consider the experimental set-up. First, the clusters formed by the laser evaporation, are cooled by the helium 
gas, and then the clusters are expanded to form a molecular beam. There is a problem in identifying the temperature 
experimentally. In the present study we assume that the clusters are in thermal equilibrium. The density of clusters 
in beam is low, so that the clusters may be assumed to be isolated beyond the equilibration zone. Therefore each 
cluster stays in a certain quantum state in the beam. Finally, the clusters enter into Stern-Gerlach magnet and are 
deflected by the interaction of the gradient of magnetic field and the magnetic polarization of super-spin induced 
by the magnetic field. At the entrance of the magnet, strength of the field changes gradually in time, and a time 
dependent interaction for the super-electron spin causes a transion of the initial quantum state to other quantum 
states. If the time dependence is sufficiently weak compared with coupling of the spin to other modes, the transition 
probability to other modes can be neglected. This is called the adiabatic approximation. In the present work, we 
calculate the profile and the magnetization with this assumption. 

If the electrons were completely decoupled from other degrees of freedom such as rotational motion, the deflection 
profile would be a flat horizontal distribution independent of the field strength. But this is not actually the case for 
the observed profiles. A small coupling of the magnetic moment to internal coordinates of the cluster gives rise to spin 
relaxation, making the profile different from the flat distribution. Hence it is important to make clear how various 
couplings produce observed deflection profiles. 

The simplest model is superparamagnetism in which the population of the magnetic states are proportional to a 
Boltzmann factoru. In other words, the cluster rotation plays a role as a heat bath for the super-electron spin in 
the magnetic field. In practical analysis for extracting the giant magnetic moment, the Langevin formula is widely 
employed. It assumes equilibrium with a thermal reservoir at a temperature which is the same as the source of the 
cluster beam. However, it predicts a rather sharp deflection profiles which is quite different from the broad profile that 
are often observed. Hence, the superparamagnetic model seems to be too simple for the analysis of the experiments. 

Another simple model is locked-spin modfiLiri which the super-electron spin is frozen to the intrinsic orientation of 
the cluster, which of course is free to rotatea E3. This model seems successful in reproducing the small peak observed 
near zero deflection angle, which experimentalists call "superparamagnetism" . But it is applicable only to Gd clusters 
and not general. Furthermore, the model ignores the angular momentum of the super-electron spin, which is known 
to give recoil effects in the Einstein-de Haas effect. 

We propose an intermediate coupling model as a method to extract the giant magnetic moment from the deflection 
profildilJ. This model covers the superparamagnetic and locked-moments models as weak and strong coupling limits, 
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respectively. This paper is organized as follows: the intermediate model is explained and formulated in Sec. 2, 
numerical results of profiles and magnetization are presented for various strength of coupling, strength of magnetic 
field and temperature in Sec. 3. The conclusions and discussions are presented in Sec. 4. 



II. INTERMEDIATE COUPLING MODEL 



In the present model it is supposed that all the electron spins are aligned in the same direction through the exchange 
interaction. Thus, the electron spins are in stretched coupling states having a giant total spin S — n s N, where n s is 
the number of spins participating the magnetic moment an atom. Since the magnetic moment is proportional to the 
total spin, the cluster has a single giant magnetic moment expressed as /i = g s S, in terms of electronic gyromagnetic 
ratio g s . 

Consider as a typical case the Feioo cluster at temperature 100K . This has a spin value of S ~ 100 and a thermal 
rotational angular momentum R ~ 600 in-,units of h. These large values of S and R justify a classical treatment of 
the problem. In fact, in a previous papeiu the classical treatment allowed the problem to be reduced a simple and 
transparent calculation. 

this was examined for a simple case, by utilizing adiabatic invariance which makes the calculation easy and trans- 
parent. 

However in a more general case, in which a fluctuation of spin orientation with respect to the cluster is significant, we 
are unable to make a simple classical treatment. Furthermore, in the classical treatment, the spin angular momentum 
is ignored, and thus angular momentum conservation is violated. Since this is a doubtful approximation when S is 
comparable to R, we follow a fully quantum mechanical treatment. 

Our Hamiltonian H is expressed as sum of three terms 



where the terms are defined as follows. The first term, H rot , stands for the rotational energy of the cluster which is 
expressed as 



where Ri's represent operators of three angular momentum components referred to the body fixed frame, and J7i's ex- 
press principal moments of inertia. The vibrational modes are not taken into account, because the Debye-temperature, 
(for instance, 500 [K] for iron) is much higher than the source temperature. In other words, the rotational motion is 
considered to work mainly as heat bath in the spin relaxation. 

The second term, Hcoupi, expresses a coupling potential between the cluster and the super-electron spin, which is 
originated from the crystal magnetic anisotropy energy caused by molecular or crystal -fields. The simplest form of 
the energy is the uniaxial magnetic anisotropy which, has already been examined in refill Here we assume that the 
clusters have an internal cubic structure and the potential has cubic symmetry. This is true for observation of the 
direction of easy magnetization being [100], [010] and [001] for iron and nickel. The anisotropy constant is measured 
in the form 



where a 2, s express the direction cosines of super-electron spin. The measured values in the bulk are K\ = 36 
[mK/atom] and K2 = 13 [mK/atom] for iron, and K\ = —4 [mK/atom] and K2 = [mK/atom] for nickel. We 
consider only the first term in the present calculation. In order to formulate the anisotropic interaction, let us start 
with the above classical picture for angular momentum variables, in which S'^s commute. 



H — H Tot + -ff COU pi + H t 



mag: 



(1) 




(2) 



E a = Ki(a\a\ + a 2 a 3 + a 3 a i) + K2a\a\a\ + ■ ■ ■ 



(3) 




= -u {S[ 4 + + S' 3 4 - 6(S[ 2 S' 2 2 + S'fS'f + s£s[ 2 ) 




(4) 
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Here (S) stands for the spherical harmonic proportional to Yx^S); $• is the spin component with respect 
to the three axes in the body-fixed frame, Si denotes the spin component referred to the laboratory system, and 
u = 2K 1 /n^N 3 . The electron spin prefers the direction of 4- fold axes of the cubic symmetry. When the direction of 
the spin is along the 4- fold axes, the value of H coup \ has minimum energy — |itS" 4 . And when the direction of the spin 
is along the eight axes of (±773 ± A= ± represented by body- fixed coordinate, the value of H coup \ turns out to be 

energy maximum y|u5' 4 . In quantum Hamiltonian, S' are operators and we are interested in their matrix elements. 
The reduced matrix element obtained through the Wigner-Eckart theorem is expressed by 




(So\\ S \\S )= - J \"' (5) 



For example, for L = 1, one obtains the familiar form ^ Sq(Sq + l)(2fio + !)■ The interaction between the rotor and 
the super-electron spin in Eq. (Q) is given by 

H coupl (s,n) = J2 [s 4 ] 4 m v 4 „: K (n)A K , (6) 



where k takes the values of only and. ±4 and A±4 — y^JliAo = — u^/32/35. The wave function of rotor can be 
expanded in terms of the 2?-functionai3, 

^(Q)=J2fRkV%(Q). (7) 

k 

Since the Hamiltonian i? C oupi is scalar, the total angular momentum I = R+ S is still a good quantum number for the 
first two terms of the Hamiltonians, H mi + H coup \ . Accordingly, we select the base labeled by total angular momentum, 
and magnetic quantum number I z . The basis is obtained by angular momentum coupling of the 2?- function T>^ ik (£l) 
and \Sa) to \IM). Therefore, the total wave function of the rotor coupled with the super-electron spin is expressed 
as 

<JW - £ (R»Sa\IM) -D* k (Wm l^> > (8) 

Rkcr 

where v represents an index specifying states having the same IM. The matrix element of the coupling term between 
the bases in Eq. (0) is estimated as 



(ffcoupl)kfl'fe' = V(2i*+l)(2S+l)(-l) s+fl - 1 

A' K (RkAn\R'k') W{RSR'S; 74), (9) 
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where ^/2Sq~+1A / ii = A k (Sq\\ S \\Sq). In our calculation we use potential strength parameter u' — A' ±4 . Now, the 

wave function (||) and the energy in the source are determined by diagonalizing H rot + -ffcoupi- Since there remains 
the (21 + l)-fold degeneracy in energy for M, we may omit M from the label for the energy E u j. 

The third term in Eq. ([lj) represents the interaction between an external magnetic field and the moments of the 
super-electron spin 

H maK = -B- n = -Bg s S z , (10) 

We choose the direction of the applied magnetic field as the axis of quantization (z-axis). This interaction breaks 
rotational symmetry for the cluster, but the magnetic quantum number M defined above is still conserved due to the 
choice of quantization axis. From (Q) and (|To|) , the matrix elements between two states are calculated as 

{9 v iM\H mae \^u'i'M) = -Bg s V2FTT{-l) 1 ' (I'MW\IM) h vIv . r , (11) 
where M independent part h u i u >j> is expressed as 

huivv = £ f&f£'VS(S + 1)(2S + 1)(-1) R - S+1 W(ISI'S; Rl). (12) 

Rk 
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For simplicity, we set the moment of inertia around the intrinsic axes to take the same value 



J\ — J2 — Jz — 3 ■ 

In this case, total wave functions are classified by an additional quantum number 717. = 0, 1, 2, 3 (k mod 4), since H mt 
is diagonal and H coup \ couples to only the states having fc-quantum numbers different by 4. This point is different 
from uniaxial anisotropy, in which fc-quantum number is strictly conserved due to the axial symmetry. The total wave 
functions in the magnetic field are expressed as 

\^m(B)) = E F *Lm(B) \*ui*>m) , (13) 

vl 

where a labels state in the same -Kk, M. 

To compare with the experiment, we calculate the average magnetization, which we regard as the ensemble average 
of (S z ). This is the sum of the product of the expectation value of S z and the occupation probability of each state. 
Now we proceed ahead with the important assumption that effect of the magnetic field is adiabatic when the cluster 
passes through Stern-Gerlach magnet. As variation of the magnetic field is slow when the cluster enters into and 
goes out from Stern-Gerlach magnet. Firstly the clusters are retained in the source region. In this region, as the 
cluster ensemble is in thermal equilibrium, the occupation probability of each quantum state is proportional to the 
Boltzmann factor exp(— -Ev/^/fceT). 

The magnetization of each state in Eq. ( |l3| ) is calculated as 

S z . a (uI)-n k M{B) = (§ a (vI)Tr k M I $z | ®ot{vI)-K h M) i (14) 

where a(vl) stands for the label of states connected adiabatically with the states {vl) defined in absence of the 
magnetic field. 

Under the adiabatic condition, any transition between energy levels does not occur even if magnetic field is applied. 
Occupation probability of each quantum state is not altered during the flight. Then the deflection profile is obtained 
by 

P(S ' B ' T) = ^) £ *( s -^KM)ex P (^) (15) 
where Z(T) represents partition function given by 

Z(T)=5>/+l)exp(=jgp*). (16) 

VlTTk 

The magnetization, the ensemble average (S z )en.av., is expressed by 

1 f —E j 

(Sz) en .av. = 777^7 2l {^<x{vI)ir k M I S z | ^(V/J^m) exp f - - ' ^ k 



(17) 



III. MAGNETIC SUSCEPTIBILITY 



In almost all the Stern-Gerlach experiments with clusters the magnetic field is so weak that the magnetization shows 
linear dependence. Hence, it is important in the analysis of the experiment to estimate the magnetic susceptibility. 
Besides, the magnetic susceptibility for the superparamagnetism and locked moment model has a well-known value; to 
examine the magnetic susceptibility will give us how intermediate coupling model describes the two extreme model. 
Therefore, we calculate the magnetic susceptibility for intermediate coupling model and discuss weak and strong 
coupling limit in the following subsections. 

For the present, let us discuss the general expression of magnetic susceptibility. Using Feynman theorem, we obtain 
the expectation value of S z as 

(* ; ) = - (d(kr) Hmag ) = ~w^) (iW ■ (18) 
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The magnetic susceptibility is expressed as 



d 



dg s B 
2 



d 2 



z 



EE 

vIM v'l' 



Bg s =0 
\(^U'I'M\S Z \^ V IM 



d{Bg a y 

2 



s=0 



E P i — E, 



v'l' 



(19) 



where Ai5 cn av . and Z means the ensemble average of energy shift and partiation function, respectively. 

We will deal with a certa in limit of susceptibility ([l9]) in the following two subs ection s. Strong coupling limit will be 
discussed in Subsect. Ill A. The weak coupling lim it wil l be treated in Subsect. [II B. Susceptibility in intermediate 



coupling will be calculated numerically in Subsect. Ill C 



A. Susceptibility for Strong Coupling 

We first consider the energy eigenvalue and the eigenstate of -ff rot + H coup i in the strong coupling limit. The 
angular momentum of the rotor R is mixed, while total angular momentum / and the z-component of total angular 
momentum, M are conserved in the absence of the magnetic field. We rewrite the Hamiltonian of the present model 
in terms of the total angular momentum, / 

1 +S \j 21 ^ + H ™&&'). (20) 

The third term of the numerator is the Coriolis term, which couples the degrees of freedom of super-electron spin to 
one of the rotor. The Coriolis term is not taken into account in the locked moment model, because the super-electron 
spin is not included as a dynamical variable. On the other hand, it must be noticed that the intermediate coupling 
model includes the Coriolis term. As will be seen this term contributes to the magnetic susceptibility. Since -ffcoupi 
gives a dominant contribution to the energy eigenvalue of H lot + H coup \ in the strong coupling, we treat the Coriolis 

term as a perturbation. The coupling Hamiltonian H coup i(S') is written in terms of the intrinsic component only 
(see Eq. (|)). Then it is advantageous to select a direct product the eigenfunction of total angular momentum and 
its z-component and eigenfunction of super-electron spin with respect to intrinsic frame; 2?jy^(f2) ® Yl^ g ak \ Sa^ as 
the basis. The energy eigenvalue of the Hamiltonian neglecting Coriolis term is 

^(1(1+1) + S(S + l))+vfE%, (21) 

where u'E^ is the energy eigenvalue of the coupling H coup \. Without the Coriolis interaction, the energy spectra is 
of rotational band with band head energy u'E^ ■ 

Let us now discuss the ground band of the present model in the strong coupling. Figure [j] illustrates the numerically 
evaluated energy eigenvalues of -ffcoupi in the space \Sa). The eigenstate of H coup i must belong to a certain irreducible 
representation of the point group O. The irreducible representations A\, E, T\ for S = 4n, E, T2 for S = 4n + 2, 
and irreducible representations T\, Ti for odd S appear in the lowest energy region. Considering the dimension of these 
representations, namely 1,1,2,3 and 3 for A\, A2, E, T± and Ti respectively, the bunch of the levels contain six states. 
The six fold and the eight fold approximately degenerate states appear at the lowest and highest energies, respectively. 
This fact implies that the direction of super-electron spin is localized to the six (eight) directions corresponding to 
the potential minima (maxima) in the lowest (the highest) bunching states for large S. Therefore, these states are 
approximated by 

R l \SS), (22) 

where R l stands for the operator corresponding to the rotation from 3rd axis to ith direction of the potential minima. 

For the present, we shall take into account the Coriolis term in first-order perturbation theory. According to 
first order perturbation theory of degenerate states, we solve the secular equation to determine the energy shift and 
perturbed states. In this case, each level in the rotational band is 6(27+ l) 2 -fold degenerate: The quantum number 
M, K have the values —I < M K < I and the ground state of super-electron spin is 6-fold degenerate. The dimension 
of the secular equation is decreased by a factor (21 + 1) because the motion under the total Hamiltonian conserves M. 
As we discussed before, the direction of super-electron spin corresponds to one of the potential minima for the ground 
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state. We choose the direction of quantization axis as one of the potential minima. This selection of quantization axis 
yields the Coriolis matrix element as 

i ■ s = iiSi + (i\s\ + As\) (23) 

Ii = R l hR l \ l\ = R i I±R i \ 
Si = R?SsR^ , S% — R i S±R^ . 



Only the first term of Eq. ( |23[) is involved in the calculation of the matrix element of Coriolis term between the parallel 
direction of super-electron spin. The second terms do not contribute the Coriolis matrix element between parallel or 
antiparallel direction but orthogonal direction. Actually, we can estimate the matrix element of raising and lowering 
operators of super-electron spin 

r 0, for 9 = 

(SS\S±R'{6)\SS) = \ cx 7h/25(l/2) s+1 , for 9 = § (24) 
I 0, for 9 = vr 

where R'{9) is the rotation through angle 9 around the axis perpendicular to the 3 axis. We find from Eq. (^4|) that 
the second term of the Coriolis term ( Eq. (|23|)) is vanishingly small for large super-electron spin. This allows us to 
approximate the Coriolis term as 

/ • S ~ hSi. (25) 

We can solve the secular equation for approximated Coriolis term, Eq. (p5|), for each degenerate rotational state. The 
wave function results in the eigenstate of P ; 



\IMN0 = \IMKi) = R\SS)R^^-V{ IK (n) (N < 6). (26) 

where £ specifies the energy levels split by the Coriolis term. For the ground band (N < 6), £ corresponds to K, and 
N is labeled by the index of the direction of super-electron spin i. The first order energy shift is obtained as 

{l'MK'j\2I ■ S\IMKi) = 2h 2 KS5 K ,K>- (27) 

We neglect the coupling between ground band and excited band because the energy difference between these bands 
are large in the strong coupling. Higher order perturbation would take into account this coupling and would give the 
deviation from the locked moment due to the finite strength of coupling. But we do not discuss the deviation in this 
paper. 

The energy eigenvalue of ( pcf ) is divided into three parts 

E IM n 6 = ^(1(1 + 1) + S(S + 1)) + u'E% + (28) 

where is the energy shift by the perturbation of the Coriolis term and is equal to Eq. ( p7| ) for the ground band 
(N<6) 

The susceptibility is evaluated through Eq. ( |19| ) as 

= -1 V V \(rMN'£\S z \IMNtj)\ 2 ( EimnA 

x ~~ z rk* z h w±i^±m _ _ u , (e a _ e a } exp r k BT ) 

(29) 

For the single rotational band with (21 + l) 2 fold degeneracy, levels of which the angular momentum are larger than 

I C (T) — J Q.2JkftT '/h 2 contributes only less than 10 percent of the partition function. We call I C {T) cut-off angular 
momentum. When the band head of excited bands are much larger than the energy of the cut-off level, that is, 

u'(E^-Et)» h2{I ; { p } \ (N>7) (30) 



G 



the higher bands do not have considerable occupation probabilities. Therefore, it is possible to describe the suscepti- 
bility with the ground band only. Substituting energy of the ground band obtained through the Eqs. (|27j) and (2S|) 
and the matrix element of S z evaluated by the wave function of the ground band Eq. (|2^) into Eq. (|9|), we get 



X 



4 E E tf- 



SwS 2 §^ i \(IMW\I'M)\ 2 \(IKW\I'K)\' 



Z ^ ^ 7(7+1) -/'(!' + 1) 



I h 2 (I(I+l)+S(S+l)+2KS) . uE A \ 

x exp ^ ^ —J . (31) 

This is different from the susceptibility for locked moment^ only in the Coriolis term. Since the average angular 



momentum is about I ~ 600, one can safely evaluate Eq. (31), treating 7, K as a continuous variable and replacing 
the sum by an integral. Thus we get 




f^dt \ - V 2 ^ (- 1 )" 2 " (acnyn+l 

dt] ~^3h* (2n + 3)!! ( ^ } 



n=0 



n 2 s 2 \ 2 \ 2h 2 s 2 



< ( 32 ) 



9 k B T \ 52Jk B T 35\2Jk B TJ "') 9 k B T ' 



where (3 = 2 j k T . The leading order is the same as the locked moment susceptibility. The higher orders mean the 
correlation by tne Coriolis term; in other words, a recoil effect due to the angular momentum of the super-electron 
spin, or Einstein-de Haas effect. We find from Eq. ([32]) that this effect reduces the magnetic susc eptib ility compared 



with one of the locked moment. This fact will be substantial in numerical calculation in Subsect. pi Q . We close this 
subsection with evaluating the condition of neglecting the recoil effect. For finite temperature, the cluster is excited by 
the thermal rotational excitation. For the rotational band with (21 + l) 2 -fold degeneracy, the most major population 

of angular momentum is I eff (T) = sj^j^. Therefore, the typical energy splitting by the Coriolis coupling becomes 

2 SI ■ ■ (T) 

h — j ' . We can neglect the recoil effect under the high temperature condition 

fc B T» JP-A (33) 



B. Susceptibility for Weak Coupling 

In this subsection we discuss the susceptibility for the weak coupling and high temperature limits. Firstly we 
discuss the energy eigenvalue and the cigenfunction of H TO t + -ffcoupi in the weak coupling. The angular momentum 
of rotor R is approximately conserved due to the weak coupling. In the absence of coupling, the energy eigenvalue 
is a single rotational band with (2R + l) 2 fold degeneracy. The coupling removes the degeneracy; the first order 
energy shift is determined by solving the secular equation in degenerate space. Therefore, the energy eigenvalue is 

described as - + AEjx(r) where A(i?) specifies states adiabatically connected with R, and AEjx(r) stands for 

the eigenvalue of the weak anisotropic coupling. The unperturbed base is expressed by 

\IMRk) = ^{RfiSallM^Sa) (34) 

The coupling is treated as the first order perturbation to degenerate states. The perturbed wavefunction is expressed 
as 

*A(fl)/M = £ (R(iSa\IM) V* k {Q)$ R) \Sa) , (35) 

ka 

where the coefficient c k is determined by the diagonalization of i7 coup i in each degenerate rotational level. 

Now we evaluate the magnetic susceptibility for weak coupling by employing the Eq. (|l9|). Equation ([l9]) reduces 
the calculation of magnetic susceptibility to the evaluation of second order energy shift. The energy denominator 
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between different R levels are much larger than the one of the same R in the weak coupling. So, we neglect the 

contribution from different R in the perturbation. The orthogonality condition X)fc c fe^* c fe = ^xy yields the 
matrix element of Zeeman interaction between the weak coupling base, Eq. (|35|), as 



*A(R)7iif|£|*v(*)i'Af ) - <Sa,a'V2I 7 TI(-1) / ' (7'M10|7M) 



xy/S(S + 1)(2S+ 1)(-1) R - S+1 W(ISI'S;R1). (36) 



Under the approximation, Eq. (19) becomes 




^X(R)Im\Sz\^X(R)PM ) i 2 E fl+ Ag j'A W \ ti z R(R+l) 

e , S/j = i— ^. (37) 



AJ VM«) AE 

Since AEj'x(R) is small because of the weak coupling, we can expand e fc B r up to first order as 1 

Rearranging the summention, we get 



I' MR) 




^X(R)Im\S z \^X'(R)I'M 



(38) 



Putting Eq. ([36j) into Eq. (pq) , and expanding for i?, Therefore, we obtain the magnetic susceptibility 

E £ =q ( 45(5 + l)(25 + l)^ + Q(J ^ exp ^±1) 



2Jk B T 



Elo ((2i? + 1) 2 (25 + 1)) exp 



(39) 



At high temperature limit many levels contribute to the summention over R. Therefore one can evaluate Eq. (|3S|), 
treating R as a continuous variable and replacing the sum by integral. Thus, we get 

9/cbJ 

ft is observed in this subsection that the locked moment susceptibility is obtained in the weak coupling and high 
temperature limit. One should note that this does not mean the super-electron spin is locked in the weak coupling. 
In fact, taking into account the axial deformation, the susceptibility of weak and strong coupling becomes different. 
In the weak coupling, the energy level is dependent on the axial deformation through the k quantum number: 

E IMRk = + !) - fc2 ) + ^jf. (41) 



The matrix element is independent of the deformation, because the Eq. (36) is independent of A. Therefore, the leading 
order of 1/ksT expansion of magnetic susceptibility is deformation independent even if the energy is dependent on 
k, and is 2S(S + l)/9 even in an axial deformed cluster. On the other hand, in strong coupling the energy level for 
the ground band expressed as 

Exmk = ^(1(1 +1)-K* + KS) + ^ 2 (42) 

If we neglect the Coriolis term, this is the same as the one of the weak coupling, apart from the degeneracy. The matrix 
element for axial deformed cluster is not different even in the strong coupling. Hence, the magnetic susceptibility of 
the axial deformed cluster is same as Eq. (|3l| ) apart from the energy. Since the matrix element divided by the energy 
difference 1(7 + 1) — /'(/' + 1) is dependent on the deformation, the magnetic susceptibility becomes deformation 
dependent. Therefore, considering the axial deformation, the susceptibility for the weak and strong coupling limits 
are different. The two limits happen to be the same value when the value of J[s are equal. This equality does not 
mean the super-electron spin is locked in the weak coupling. 
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C. Intermediate Coupling 



We numerically calculate the magnetic susceptibility for intermediate coupling in this subsection. In numerical 
calculations of the deflection profile and the magnetization, a main task is diagonarization of Hamiltonian matrices of 
large dimensions. In the source area, the partial Hamiltonian -ff mag is zero, and therefore the total angular momentum 
is conserved. The dimension of matrices to be diagonarized is approximately S x I C ^{T). The dimension of matrices to 
be diagonarized is of the order of 10 5 . In the magnetic field where states having 1/2 different are mixed, the number 
of dimension is magnified by 1/2, and becomes eventually of the order of 10 7 . It may not be feasible in numerical 
calculations. Let us use a similarity transformation in order to scale down the angular momenta. In Eq. (|^), the 
matrix element is approximated as 

(H ) -A' V( 2 ^+l)( 2 ^ + l) ,4 (f) W 4 (f) x u ^ 
{n couv \) R K,R'K' — A K — ; «K/i(W<Vo(y2,) 



with 



2 

and two angles, 9\ and 62, are defined respectively as 



2R+1 



R = and A* = Rl - R, (44) 



K . , — K + K' 

cos 9\ = \ = with K = , 

W R 2 



and 



O2 (R-S + I)(S -R + I) 

cos— = 1/ = . (45) 

2 V 4RS ' 



These angles are invariant under a scale transformation defined as 



' r,r' r_\ /R,R'R_\ 

k, k' k = 77 1 K, K' K . (46) 



\s,i,m/ \ S,I,M J 
In a similar manner, the matrix elements given in Eq. (fTD) are expressed approximately as 



X^vkjcIM I -ffmag | ^v'I'tt k m) 



21 

Y.^KfRK{-) I '- I dl I - I >{e 5 )d 1 I _ II0 {6 i ) (47) 

RK 



with 



cost's = \ — 



cos-^/ — S + R)_(S-I + R) 

The two additional angles, #3 and 64 are also invariant under the scale transformation. Eventually these four angles 
are invariant with respect to the scale transformation in Eq. ([46|). The four d- functions are all smooth functions of 
the four angles. 

In the numerical calculation we scale down the super-electron spin quantum number S. To make S 2 invariant under 
the similarity transformation, Eq. (^|), the Plank constant is changed according to rj. In other words, we can scale 
down the angular momentum quantum number by adjusting Plank constant as 



The temperature fee? 1 , coupling it', and magnetic field parameter g a BhS, which have dimension of energy, are invariant 
under the similarity transformation. In the following, the unit of these energy parameter is selected as ^^j- which is 
also invariant under the similarity transformation. 

In the numerical calculation, we set the cut off of super-spin s — S c and the cut off of total angular momentum 
/„. T) is decided through the S c and Eq. ( pfcf ) as 77 = S c /S. 

We are now ready to calculate numerically the susceptibility for intermediate coupling. It is already known by the 
analysis in the above two subsections that the susceptibility becomes the locked moment in the limits. To achieve 
locked moment limit, the ground state of the H conp i has to be approximately six fold degenerate. When we diagonalize 
-ffcoupi in S c = 10 (see Fig. [l]) three energy levels having 1-fold, 3-fold, and 2-fold degeneracy from lower energy to 
higher make a bunch around the ground state. The energy differences of these states are about 0.01. In contrast to 
this, we find from Fig. [l] that the energy difference between bunch of the ground state and the first excited bunch 
is about 1.0. The energy spacing between six states around the ground states are much smaller than the energy 
spacing between the ground bunch and the first excited bunch. Therefore, it seems reasonable that the ground state 
is approximately 6-fold degenerate for S c = 10. 

A truncation of the angular momentum is necessitated due to computational limitations. We select I max = 80 in 
the calculation. As mentioned previously, we do not confident about the calculation of partition function out of the 
reliable range of temperature. Hence, the numerical calculation of partition function is reliable only I C {T) < J ma x- 
This restricts the range of temperature as k-gT < 20.6. We choose 20 for the temperature. 

We display the magnetic susceptibility of our calculation in Fig. ^|. The first point to be discussed is whether our 
calculation reaches the two distinct limits discussed above. If the coupling is much smaller than the temperature, the 
susceptibility decreas es an d converges to 2/9. Therefore our calculation agrees with the analysis for weak coupling 



discussed in Subsect. IIIB 



Let us discuss the strong coupling limit. In Subsect. Ill A we discuss the condition of temperature and coupling in 

re condition for the coupling and for the temperature are 



which the susceptibility attains the locked moment value. T 

expressed in Eq. (|3Cl) and Eq. (|33|), respectively. These two conditions for S c = 10 become 



u 

>3.1, and fc B T>4. (50) 



Therefore, we can regard ~ 9 for k^T = 20 as the strong coupling region. However, in Fig. ||, the susceptibility is 
not 2/9 even if the coupling becomes larger than 9. In particular, one might not readily believe that the susceptibility 
is smaller than 2/9. One should recall that we evaluate the strong coupling limit of susceptibility including the effect 
of Coriolis coupling as Eq. (|32|). We find from this expression that the susceptibility is decreased by the recoil effect. 
We calculate the susceptibility for k-^T — 20 through Eq. (p2[). The result is illustrated in Fig. ^ and is close to 
the numerically computed susceptibility. Nevertheless, the numerically computed susceptibility shows some deviation 
from the value of Eq. (|32"|). A possible reason for this is that the quantum fluctuation between different potential 
minima of the coupling decreases the susceptibility. 

We find a peak at = 0.55 in Fig. ||. This peak is accounted that the energy splitting of anisotropic interaction 

2u' is comparable to the energy splitting by Coriolis term jj4I c ffS. The width of energy splitting of anisotropic 
interaction corresponds to the energy difference between potential maxima and minima i.e. In practice the width 
becomes smaller due to the quantum fluctuation and is found 2u' for S = 10 from Fig. [j]. Then the peak is expected 
to be, 

2u '-S 4/ ^ ~ ^ 2 V^^ - 45 - (51) 

In fact, this condition is consistent with the observed position of peak. 

Let us move to the next point. As mentioned in introduction, one often assumes the superparamagnetism, in which 
susceptibility is ^S(S + 1), for the clusters to analyze the experiments. From Fig. || that the superparamagnetic limit 
is not reach in any range of coupling strength. 

The last point is the temperature dependence of the magnetization. It is reasonable to consider that the susceptibil- 
ity decreases as the temperature increased because of thermal fluctuations. A quite opposite dependence is reported 
in Refill. However, our calculations does not show such a behavior in any range of coupling. 
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IV. PROFILES AND MAGNETIZATION 



As mentioned in the previous section, a truncation of the angular momentum is necessary due to the computational 
limitations. The procedure of calculating profile includes the large dimensional diagonalzation of Hamiltonian H mag . 
To make it practicable, we choose a truncated value of J ma x = 26, and the magnitude of super-electron spin S c = 10. 
As discussed in section III, this choice narrows the reliable range of temperature. The reliable range of temperature 
for this cut-off becomes k^T < 2. Since high temperature is one of the conditions for the strong coupling limit, we 
should select as large temperature as possible. Eventually, we choose k^T — 2 which is maximum for the reliable 
range of 7 max = 26. 

Firstly we pay attention to the magnetic field dependence of the deflection profile. Figure || shows the deflection 
profile for three different strengths of the magnetic field g s BTiS = 0.5, 2, 10, respectively. Many experiments show that 
peak deflection is always to the strong field. Our calculation also exhibits deflection of a similar nature. By solving 
the classical equation of motion, De Heer et al. have demonstrated in ReJld that the coupling between super-spin 
and cluster body causes the deflection behavior. We will have a closer examination of this statement in terms of our 
quantal model. 

In order to make our discussion transparent, we discuss the weak coupling limit first. We can see from Fig. |3| 
that the profile of the weak coupling is a flat but slightly sloping distribution in the weak magnetic field. Figure 
[|, which shows the energy levels of the weak coupling v! = 4.47 x 10~ 3 for 5 C = 2,M = 0,7Tfe = and R ~ 3 
or 4, helps us to understand the distribution in terms of the energy level. In the absence of magnetic field, the 
weak coupling removes 15- fold degeneracy of unperturbed rotational levels; Energy levels result in a bunch of states 
around the unperturbed rotational levels. The applied magnetic field further splits these levels, and rearranges 
each one more likely to be the eigenfunction of S Zl when the applied magnetic field becomes stronger than the 
coupling. In other words, the super-spin precesses about the direction of magnetic field independently from the 
cluster (decoupling) for g s BTiS > u' . Although the magnetization increases proportional to the magnetic field before 
the decoupling, occupation probabilities of the same R levels become almost equal because of the weak coupling. Then, 
the magnetization remains small and the profile becomes flat, rectangle like distribution. When the Zeeman splitting 
equals to the typical energy difference of the rotational levels: 2g s BhS — 2Ti 2 R cS /2J, the pseudo-crossing between 
different rotational levels occurs. To put it another way, the pseudo-crossing occurs where the Larmor precession 
frequency c<jl is comparable to the cluster rotation frequency w rot ; wl — w ro t ■ That is to say, the pseudo-crossing leads 
to exchange the occupation probabilities between those levels. Therefore, the magnetization increases and profile 
develops a peak at (S z )/S = 1 like the Boltzmann distribution. 

The above discussion leads us to divide the mechanism of magnetization into two types, that is, the magnetization 
by the process of decoupling and pseudo-crossing. When the coupling is weak, these two are separated definitely. 
This makes the magnetic field dependence of magnetization peculiar as shown in figure The magnetization linearly 
increases as the magnetic field in the process to decoupling; then, remains steady by the decoupling of super-spin; 
finally, increases suddenly because of the pseudo-crossing. Since the magnetic field region of the decoupling and 
pseudo-crossing can not be divided into the intermediate and strong coupling, the peculiar behavior of magnetization 
disappears. 

Let us now turn to the strong coupling. The calculated susceptibility (Fig |J) indicates that v! — 11.2 is in the strong 
coupling region. However, the calculated profile for v! = 11.2 in figure ^ is not identical with the locked moment 
profile. We attribute this to the Coriolis term, which is important at the temperature of our ensemble, k-gT = 2. A 
higher temperature was possible in the calculation in refclll because simpler anisotropy term there permitted smaller 
dimension Hamiltonian matrices. 

Finally, let us move to the discussion of the profiles in intermediate coupling. We characterize the profiles by the 
cumulants. Since the cumulant higher than the second order vanishes for the Gaussian probability distribution, the 
characterization of the profiles by cumulant clearly shows us the deviation from Gaussian profile. In our calculation 
the third and the fourth order cumulants are much smaller than the first and the second one, so the Gaussian profile 
gives a reasonable description. 

We now discuss the first (Fig. |?j) and the second order cumulant (Fig. |^) in more detail. These correspond to 
ensemble average of the magnetization and width of the profile, respectively. We can estimate the second order 
cumulant of the locked moment profile in the absence of magnetic field and the rectangle profile as | and S< '^2 1 ' > , 
respectively. The general trend of the first order cumulant looks like more or less one of either the superparamagnetic 
or locked moment. Although, we introduce the present model for describing the intermediate behavior between 
superparamagnetic and locked moment, the magnetization of intermediate coupling is smaller than one of the locked 
moment. This is because the Coriolis term decreases the susceptibility, as we discussed in Sec. III. We may note, in 
passing, that the anomalous behavior in the weak coupling, discussed earlier, should appear when the energy splitting 
by the magnetic field comparable to typical energy difference of rotational level, that is, 
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2g sB hS^^^ * g s BhS^U h2S y T -0.14 (52) 
y 2J y S]j 2 J y ' 

But it is not seen in the figure 0. 

The second order cumulant is sensitive to the coupling. The stronger the coupling becomes, the narrower the 
profile. The cut off angular momentum in this calculation is rather small than the perturbation. We calculate 
cumulant for higher temperature by applying the perturbation technique. Figure ^| shows the second order cumulant 
for /max = 80, k#T — 20 as the function of the coupling strength. 

The second order cumulant becomes closer to the strong coupling limit S 2 /9 around the coupling u' — 3.5. It 
decreases almost linearly in the region of the coupling beyond 20 due to the tunneling between different direction of 
super-electron spin. Since the perturbation theory is valid for the magnetic field weaker than the coupling strength, 
the super-spin do not decouple from the cluster. Therefore, the profile is not the rectangle profile. We can estimate 
the second order cumulant of weak coupling limit. We find the magnetization of each level in this limit from the 
diagonal part of Eq. ; 



(S z ) = \/27TI(-l) / (IMW\IM) y/S(S + 1)(2S + l)(-l) R - s+1 W(ISIS; Rl). (53) 
The second order cumulant for weak coupling is estimated as 



RIMX(R) 



(54) 



We treat R as a continuous variable as in the Subsect. [II A| and get 



(•^z) cumulant — g^(^ + -0 (^5) 

Actually, in Fig. ^, the second order cumulant is not close to the one of the rectangle profile 1/3, but close to 
5c( 9 % +1) ~ 0.122. 

To make the profile rectangle, the super-spin is decoupled from the cluster because of the strong magnetic field 
compare with the coupling; 

g s BHS > u'. (56) 

Furthermore, pseudo-crossing do not occur at least for the levels having a typical angular momentum. Because, 
pseudo-crossing changes occupation probability of the levels. 

h 2 

g s BfiS < ^j2R cff (57) 

However, in the perturbation for the magnetic field, the magnetic field is always smaller than the coupling. Hence, 
we can not describe the rectangle profile by perturbation for the coupling. 

To describe the regime around the rectangle profile, we calculate magnetic field or coupling dependence of (S z ) 2 
by perturbation of the coupling under the two assumptions discussed above. The detailed account of calculation is 
presented in appendix. Here we give just the result Eq. ( |A16 ) . 



(s z ) 2 = ls(s + v s ~ ( 



3 v 1 15 x 9 2 \Bg s h 

190 + 395(5^ 



2 



9x15V 693 J -i^ 5(5 + 1) jU^J ' (58) 

The higer order cumulant shows the deviation from the gaussian profile. As a detail of the calculated higher order 
cumulants from obtained profiles are not shown, the result is much smaller than the second order. When one analizes 
the experiment in more detail, the higher order cumulants may bring about important imformation. 

For example, the third order cumulant (Fig. ^ characterizes asymmetry of the profile. In the whole range of 
magnetic field, the third order cumulant is smaller than the second order and is smaller as the coupling becomes 
stronger. The profiles are always symmetric with respect to (S z ) = in the absence of magnetic field. When the 
magnetic field is applied, the asymmetry grows on account of the time reversal symmetry breaking. The third order 
cumulant have a peak at g s B/k^T ~ 2. In the strong magnetic field, the profiles have a narrow peak at = 1 
like Fig. |^; The profiles are more likely to be symmetric but still not quite. In other words, the third order cumulant 
decreases but is not equal to zero. 
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V. CONCLUSIONS AND DISCUSSION 



Using an intermediate coupling model, we have studied the magnetization of ferromagnetic clusters in a Stern- 
Gerlach magnet. In this model the super-electron spin is free but couples to the cluster ions through an anisotropic 
potential. This model is expected to describe the intermediate behavior between the superparamegnetic and the 
locked moment. In evaluating the profiles or the magnetization, we assume that the variation of the magnetic field in 
entering the magnet is slow in time, i.e. adiabatic. Hence, any transition between the quantum states is suppressed; 
the occupation probability of each quantum state is determined in the source area where the magnetic field is absent. 

We examined the magnetic susceptibility of the present model applying perturbation theory. Especially, the mag- 
netic susceptibility in the strong and in weak coupling limits are discussed analytically. We expected the intermediate 
coupling model to approach the locked moment behavior in the strong coupling limit. However, there is a crucial 
difference between the strong coupling limit of the present model and the locked moment model. In the present model 
the super-electron spin degree of freedom is treated explicitly and the Coriolis term arises due to the conservation 
of total angular momentum, while in the locked moment model the Coriolis term is neglected. Consequently, the 
magnetic susceptibility is different from, and smaller than the locked moment. In the weak coupling limit, we expect 
the susceptibility of the present model to become superparamagnetic value. But the susceptibility in the present 
model obtained by the perturbation theory is not the superparamagnetic value in any range of the coupling. 

We find that the irregular magnetic field response noticed in re£2l also exists in the weak coupling region of the 
present model. The magnetization linearly increases before the decoupling; and is saturates after the decoupling. 
Then, once the pseudo-crossing between different R levels occurs, the magnetization increases again. Before the 
decoupling, the susceptibility arises as a perturbation of the magnetic field. But after the decoupling, or once the 
pseudo-crossing takes place, the susceptibility becomes nonperturbative. Instead, we need to diagonalize the total 
Hamiltonian, which is difficult for today's computers for high temperatures and large cut off angular momentum. In 
re£3, one of the authers (GB) discussed the susceptibility after the decoupling using uniaxial coupling. They found 
that the susceptibility reaches the superparamagnetic value. Therefore, if we would calculate the susceptibility after 
the decoupling, the susceptibility would be superparamagnetic. n 

Although an anomalous temperature dependence is reported in rem, the calculated susceptibility is always positive. 
We do not reproduce the anomalous temperature dependence. 

The " superparamagnetic peak " which is seen in the profile obtained by the locked moment model is not seen in 
the present calculation even in the strong coupling limit. As discussed before, the strong coupling limit of our model 
is different from the locked moment model due the Coriolis term. Hence, the effect of Coriolis term destroys the 
superparamagnetic peak. If we would be able to calculate in such high temperature that the Coriolis term can be 
negligible, the peak would appear in the deflection profile. 

Finally, we summarize our method of analysis of the Stern-Gerlach deflection function. We calculated the cumulant 
of the profiles up to 3rd order to characterize the profiles. We found from figure that the first and second order 
cumulants are dominant. In particular, the evaluation of the susceptibility and second order cumulant by the pertur- 
bation technique gave us the analytical expression of the second order cumulant and magnetic susceptibility for high 
temperature and strong or weak coupling. One can extract the magnetic moment and the coupling strength by fitting 
observed magnetic susceptibility and the second order cumulant into the calculated second order cumulant and the 
magnetic susceptibility. 
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APPENDIX A: THE SECOND ORDER CUMULANT FOR THE DECOUPLING REGIME 

In this appendix we derive Eq. ([58]) under the assumption Eq. (|56| ) and Eq. (^?j). The ensemble average of the 
second order cumulant determined through (S z ) 2 of each level and the energy in the absence of magnetic field. Let us 
start the calculation of second order cumulant with evaluating (S z ) 2 . The assumption Eq. ( |56"| ) assures us of applying 
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perturbation theory of the coupling. Then we treat H coup \ as perturbative and H mt + ff mag as nonperturbative. The 
degeneracy remains for intrinsic quantum number k when we select 2?^ fe (f2) <8> \Scr) as unperturbed bases. According 
to perturbation theory for degenerate case, the unperturbed bases are obtained as eigenstates of the Hamiltonian 
^coupi of the degenerate space; 



(0) 

RMcrv 



2R+1 



8tt 2 



c£2>£ fc (Jl) ® jStr), M = fi + a 



(Al) 



The effect of the coupling is included by first- and second-order perturbation theory. In the perturbed wave function, 
there exists mixing with different rotational states. Now we should recall the assumption Eq. ( |57| ) in which the energy 
difference between different rotational states are much larger than the magnetic field. This justifies that we can neglect 
the mixing with the different rotational states. Therefore, the perturbed wave function is evaluated as 



(2) 

RMcrv 



u' 2 \(h conpl )X' 



-y 

2^, {Bg s h{a-a')Y 



(0) 



RcruM 



u (^ , coupl) I/(T ' |I / - ^(o) 



cr^cr' 



Bg s h((T — a 1 ) 



Rcr'vM ' 



M- 



(A2) 



with 



(0) 



u (^coupl) v i cr' ,1/cr — ' Rcr 1 V ^ A/ I coupl I ^ RcryM 1 



I* 



(o) 



u'5 vy A vR (Rij/4m\ Rfx) (Saim \Sa'), 



(A3) 



where A v r represents the eigenvalue of hf.,k> = E K ~X^i {RkAK,\R'k') and u' = A±± ^^s+f ^ ' Using the perturbed 
wave function Eq. (A2), we calculate (S z ) 2 up to second order for the coupling, 

.,,'2 



{Bg s H) 2 
where C'2 = 



C 2 , 



I (^COUpOj/'^'^cr I 

P - a) 



(A4) 
(A5) 



Let us move on the evaluation of the energy in the absence of magnetic field. We should evaluate the energy 
which is adiabatic ally connected with state in the finite magnetic field. If the coupling is absent, the unperturbed 
wavefunction Eq. (Al) adiabatically connects between states in the magnetic field and one in the absence of magnetic 
field. Hence, we start from unperturbed wavefunction Eq. (Al), then introduce the coupling up to the second-order. 

According to the perturbation theory, the first order energy shift u' 'Cyo is evaluated the expectation value of -ff C oupi 
for unperturbed state; 



/ e W = {Rfj,40\Rfj){SaA0\Sa)u'A vR . 



(A6) 



In the second order energy shift, different rotational states mix. However, as is mentioned above, the assumption Eq. 
([57]) enables us to neglect the mixing of different rotational states. Ther efore we negle ct second order energy shift. 

We are now able to calculate ensemble average of (S z ) 2 using Eq. (A5) and Eq. (A6). Expanding not only (S z ) but 
also the Boltzmann factor up to second order, we get the ensemble average of {S z ) z as 



(S z 
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k B T 
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" ^uR 



E 



RucrM 



.(1) 
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%CD a 

t vR 



exp 



exp 



-E R 
k B T 



(A7) 



E„ e uR = 0: E„ a2 ^vR — 0' because Tr A vR = K (Rk40\Rk) = 0. Expanding Eq. ( JA7| ) up to second order of vf, we 
get 
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cn.av. 



5(a 2 ) , S(C 2 ) ( v! \ 2 1 ( 5(a 2 e« ) 5(e« )5(a 2 ) 

2 



5(1) 5(1) \Bg s h 



5(1) 



5(1)2 



where 5(x) is defined as 



S{x)= ^ xexp(-E R /k B T) . 



RucrM 



Let us estimate S(x)'s appearing in Eq. (A8). 

5(1) = ^2(2S + l)(2R + l) 2 exp 



-E R 



S(cr 2 ) = J2 l S ( S + !)( 2 ^ + !)( 2i? + !) 2 ex P 



(A8) 
(A9) 

(A10) 
(All) 



Before calculating the other 5(x)'s, we evaluate J2u \^vr\ 2 - Since the trace of a matrix is constant under the unitary 
transformation, we get, 



vR 



2 = TriU^H^UU^H^U) = Tr(tf 2 oupl ) 



V(^ | (Rk\A0Rk}\ 2 + | (Rk\MRk + 4) | 2 + | (Rk\A - 4Rk - 4) | 2 ) 



fc 

8(2i?+ 1) 
15 



Substituting Eq. (A12), we evaluate the other S(x)'s 

Q( W\ - \" 8(2^+1)2^+125 + 1 / -En 



■ exp 



knT 



a , 2 (i)2 ^8(2#+l) 2R + 1 (1 + 25)(-190 + 395(5 + 1)) f-E R \ 
S{a €vr) ~^ 15 9 693 eXP {l^T 

, V- 8(2R+l)2R +1-8(25+1) 
^ 2 )-L^ g g ex 



(A12) 

(A13) 
(A14) 
(A15) 



Finally, putting Eqs. ( |AlOj ),( [ATl| ), ( |Al3| ),( |AT4l ),( [A15| ) into Eq. ( |A8[ ), we obtain the ensemble average of (5 Z ) 2 as a 
function of the coupling strength. 



(S z ) 2 = -S(S + 1) — , , 

V 1 3 V ' 15 x 9 2 \Bg s h 

4 /-190 + 395(5 + 1) 



9 x 15 
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(A16) 
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Figure 2 



0.26 



0.25 



0.24 



0.23 - 



0.22 



0.21 



-A 


















+ 

+ 4 


+ 












k B T=20 




+ 

+ 


+ 

\ 














- 


+ 


+4 








Locked n 


loment 






+ 




































Loci 


:ed mome 


<nt with C 


!oriolis te 


rm 


H- + 4f +f + 



























0.2 







2 3 4 5 6 7 

Coupling Strength u'/k B T 



8 



FIG. 2. The magnetic susceptibility calculated by the perturbation theory of the intermediate coupling model as a function 
of the ratio of coupling strength to temperature. 
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Figure 3a 
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Figure 3b 
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FIG. 5. An example of peculiar behavior of magnetization for S c 
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Figure 6 
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FIG. 6. The second order cumulant as a function of the ratio of coupling strength and temperature. The temperature is 
fixed at 20. 
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Figure 7a 
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Figure 7b 
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FIG. 7. The first order cumulant of each profile as a function of the ratio of magnetic field. The temperature is fixed at 
k B T = l and for k B T = 2 for Fig. 7a and Fig. 7b respectively. 



23 



Figure 8a 
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FIG. 8. The second order cumulant of each profile as a function of the ratio of magnetic field. The temperature is fixed at 
k B T = 1 and for k B T = 2 for Fig. 8a and Fig. 8b respectively. 



24 



Figure 9a 
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Figure 9b 
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FIG. 9. The third order cumulant of each profile as a function of the ratio of magnetic field. The temperature is fixed at 
k B T — 1 and for k B T = 2 for Fig. 9a and Fig. 9b respectively. 
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